Dynamics of entanglement via propagating microwave photons 



C. Sabi'n, 1 J. J. Garcia- Ripoll, 1 E. Solano, 2,3 and J. Leon 1 

1 Instituto de Fisica Fundamental, CSIC, Serrano 113-B, 28006 Madrid, Spai^^ 
Departarnento de Quimica Fisica, Universidad del Pais Vasco - Euskal Herriko Unibertsitatea, Apdo. 644, 48080 Bilbao, Spain 
3 IKERBASQUE, Basque Foundation for Science, Alameda Urquijo 36, 48011 Bilbao, Spain 

(Dated: May 3, 2010) 

We propose a simple circuit quantum electrodynamics (QED) experiment to test the generation 
of entanglement between two superconducting qubits. Instead of the usual cavity QED picture, we 
study qubits which are coupled to an open transmission line and get entangled by the exchange 
of propagating photons. We compute their dynamics using a full quantum field theory beyond the 
rotating-wave approximation and explore a variety of regimes which go from a weak coupling to the 
recently introduced ultrastrong coupling regime. Due to the existence of single photons traveling 
along the line with finite speed, our theory shows a light cone dividing the spacetime in two different 
regions. In one region, entanglement may only arise due to correlated vacuum fluctuations, while 
in the other the contribution from exchanged photons shows up. 

PACS numbers: 03.67.Bg, 03.65.Ud, 85.25.-j 



I. INTRODUCTION 

Quantum mechanics does not allow us in general to 
consider two arbitrary distant systems as separate^. In 
some cases there exist quantum correlations that cannot 
be generated by local operations and classical communi- 
cation between remote systems. Time enters this picture 
through two different questions. The first one is related 
to the speed bound of a hypothetical superluminal influ- 
ence which could explain all quantum correlations, esti- 
mated to be 10 4 c in a recent experiment!^. The second 
question is of a more practical nature inside the quan- 
tum theory^HIl what is the speed at which two distant 
systems become entangled? 

Quantum field theory (QFT) fulfills the principle of 
microscopic causality by which two space-like separated 
events cannot influence each othepSl and thus cannot 
be used to transfer information^^. We may then ask 
whether microcausality also sets a limit on the speed at 
which entanglement can be created between two separate 
systems. More precisely, can two subsystems, supported 
at regions (x, t) and (x', t'), become entangled while they 
are still space-like separated? Or in simple terms, can fi- 
nite quantum correlations develop before signals arrive? 
The answer to this far reaching question is yes, it is pos- 
sible. After all, Feynman propagators are finite beyond 
the light cone and even before photon arrival there exist 
correlations between the vacuum fluctuations at any two 
space-like separated events. 

In this work we demonstrate that circuit QED is ar- 
guably one of the most suitable fields to study the dynam- 
ics of entanglement between distant systems. One rea- 
son is the existence of various choices of high quality su- 
perconducting qubits, the so-called artificial atoms^"^. 
Another reason is the possibility of coupling those qubits 
strongly with traveling photons using microwave guides 
and cavitief^H^. Furthermore, those co upling strengths 
can reach the ultrastrong coupling regime^HSS^ where the 
qubit-photon interaction approaches the energies of the 



qubit and photons. In this case, the rotating-wave ap- 
proximation (RWA) breaks down and a different physical 
structure emerges. Such regimes can be activated and 
deactivated^, facilitating the creation of a fairly large 
amount of entanglement in a time-dependent way, as we 
will see in this work. 
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FIG. 1: (a) Qubits that interact via traveling photons with 
finite velocity v can be space-like (I, white) or time-like (II, 
shaded) separated, depending on the value of £ = vt/r. While 
only in II they are causally connected, entanglement may ap- 
pear already in region I. (b) A possible implementation of 
these ideas consists of flux qubits ultrastrongly coupled to a 
common transmission line, (c) With a slight modification, the 
coupling of the qubits to the line can be dynamically tuned 
via fast magnetic fluxes, $ (Color online). 
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We will discuss some of the preceding questions in the 
framework of a precise circuit QED setup, see Fig. [I] con- 
sisting on two well separated superconducting qubits cou- 
pled ultrastrongly to an open transmission line [Fig. |T]b] . 
The waveguide provides a continuum of microwave pho- 
tons propagating with uniform velocity, v, mediating an 
interaction between the qubits. Given an intitial sep- 
arable state in which only qubit A is excited, we have 
studied the evolution of correlations and related it to 
the propagation of photons between qubits. The main 
results are: (i) Outside the light cone, that is in re- 
gion I of Fig. [iji where £ = vt/r < 1, the excitation 
probability of qubit B is independent of the distance r 
to qubit A. (ii) Still in region I, entanglement between 
the qubits always takes a finite value and grows with 
time, (iii) Once the qubits are time-like separated, that 
is as soon as we cross into region II, entanglement grows 
faster than the excitation probability of qubit B and 
takes sizeable values. Result (i) is a manifestation of the 
fact that our Quantum Field Theory (QFT) model satis- 
fies microscopic causality, which formally translates into 
the vanishing of commutators associated with observ- 
ables at space-like separations, [<2(x, t), Q'(x', t')] = 
for |x — x'| 2 — c 2 (t — t') 2 > 0. Furthermore, it shows 
that two qubits which are space-like separated cannot be 
used to communicate superluminal information. Result 
(ii), on the other hand, reveals the fact that correlations 
between vacuum fluctuations at separate points can be 
established at arbitrarily short times, even though they 
are non-signalling and cannot transmit information. 

It is important to remark that the previous questions 
have been posed theoretically using model detectors^, 
two-level atomJ^', scalar fields^l and photoniPH^, yet 
no experimental test has been accomplished. However 
in this work, we show that the access to the ultrastrong 
couplings in circuit QED allows us to explore these ideas 
with very advantageous parameter ranges. 



II. SUPERCONDUCTING QUBITS COUPLED 
TO A QUANTUM FIELD 

Our setup consists of two qubits, A and B, interacting 
via a quantum electromagnetic field. The qubits have 
two stationary states |e) and \g) separated by an energy 
Ml and interact with a one-dimensional field, which prop- 
agates along the line connecting them, 

V(x) = J dk^/NuT k [e tkx a k + H.c] . (1) 

This field is described by a continuum of Fock operators 
[a/-,ajj] = Skk'i an d a linear spectrum, ui^ = v\k\, where 
v is the propagation velocity of the field and plays the 
role of the speed of light. The normalization and the 
speed of photons depend on the microscopic details. In 
particular v = 1/%/cI, where c and I are the capacitance 
and inductance per unit length. 



We consider qubits that are much smaller than the 
relevant wavelengths, A — v/Q,, and lay well separated. 
Under these conditions we can split the Hamiltonian, 
H = Ho + Hj, into a free part for the qubits and the 
field 

H Q = ^hn(cr z A + a z B )+^2oj(k)ala k , (2) 

and a point-like interaction between them 

H r = J2 d a V(x a ). (3) 

a=A,B 

Here xa and Xb are the fixed positions of the atoms, and 
d n = d x is equivalent to the dipole moment in the 
case of atoms interacting with the electromagnetic field. 

In what follows we choose the initial state \ip{0)) = 
\eg) <8> |0), where only qubit A has been excited, while 
both B and the field remain in their ground and vacuum 
states, respectively. In the interaction picture given by 
the "free" Hamiltonian Hq, the system evolves during a 
lapse of time t into the state 

|^(t)) = T[e- % ti dt ' H '^l h ] \eg) ® |0) , (4) 

T being the time ordering operator. Up to second order 
in perturbation theory the final state can be written as 

\i/>(t)) = [(l + A)|e S )+X| fl e)]®|0> + 
aj A \gg)+V B \ee))®\\) 
+(F\eg)+G\ge))®\2)+0(d 3 ). (5) 

The coefficients for the vacuum, single-photon, and two- 
photon states, are computed using the action (a — A, B) 

Si = - % lf Q (*«\d°Z \g a ) V(x a , t')dt' = - (S- )t (6) 

among different photon number states \n) , n = 0, 1, 2 . . ., 
being \n) (n\ = ^ J dki.... J dk n \ki...k n ) {ki...k n \ and 
\k) = at |0). Only one term corresponds to interaction 

X = (0\T(S+S^)\0). (7) 

This includes photon exchange only inside the light cone, 
vt > r, and vacuum fluctuations for all values of t and r, 
being r = xb — xa the distance between the qubits. The 
remaining terms are 

A= 1 -(0\T(S+S A +S B S+)\0) (8) 
C^ = (l|5^|Q),Vb = (l|5+|Q) J 
F = i <2| T(S+S A + S B S+) |0), G = (2| T(S+S A ) |0) . 

Here, A describes intra-qubit radiative corrections, while 
Ua,Vb,F and G correspond to single-photon emission 
events by one or more qubits. 
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Note that virtual terms like V, F and G, which do not 
conserve energy, are relevant only at very short times 
and are always neglected in the RWA. Here, we are in- 
terested in the short-ti me beha vior, and therefore all the 
terms must be included 9 * 22 * 23 !. Furthermore, only when 
the non-RWA terms are included, can it be said properly 
that the probability of excitation of qubit B is completely 
independent of qubit A when r > vt, which is precis ely 
the condition for microcausality discussed bef or #21221. 

The coefficients in Eq. ^ can be computed analyti- 
cally^ as a function of two dimensionless parameters, £ 
and K. The first one, £ = vt/r, was introduced before 
and it distinguishes the two different spacetime regions 
[Fig.JT^i] , before and after photons can be exchanged. The 
second parameter is a dimensionless coupling strength 
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Note that the qubit-line coupling g = d^NQ/K cor- 
responds to the qubit-cavity coupling that appears by 
taking the same transmission line and cuttin g to h ave a 
length L = A thus creating a resonator RefsP^ESl This 
formulation has the advantage of being valid both for in- 
ductive and capacitive coupling, the details being hidden 
in the actual expressions for d and N. 

Tracing over the states of the field, we arrive at the 
following reduced density matrix 



/ Pn 



PAB = ~ 
r 



pi4 



, _., , .... 
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P22 P23 
P* 23 P33 



(10) 



representing the two-qubit state in the basis formed by 
\ee) , \eg) , \ge) , and \gg) . The coefficients with the lead- 
ing order of neglected contributions are 



Pn 

P33 

pu 

P23 



\v\l 
\*\ 2 

U* A V B +0{d 4 ) 
X* +0(d 4 ), 



0{d 4 ), P 22 = I + 2Re( A) + 0(d 4 



\G\ 2 + 0(d 6 ), p 44 = 
(0|5+5+| 



\U\ 2 A + 0(d 4 ) 
0)+O(d 4 ) (11) 



and the state is normalized c = pa. 

Let us now remark the validity of the perturbative 
methods applied in this work. The leading corrections 
to C(pab) (see Eq. (12) below) come from the lead- 
ing order corrections to P23jPi1jP44 

(Eq. ([II])). In the 

case of p 23 we have p2 3 (d 4 ) = (1 + A)X* + FG* and 
P23{d 6 ) = P23(d 4 ) + Xi + X 2 , where X x comes from the 
interference of one and two photon exchange amplitudes 
and X2 comes from the probability amplitude of three 
photon exchange. A rough upper bound for these two 
terms is given by 2 |A| 3 . For pn and P44 they involve 
a number of photon emissions and reabsorptions by the 
same atom or by the other, giving a term pnp44(<i 6 ) = 
|J7a| 2 |Vb| 2 + A\ + A 2 , where rough upper bounds to A\ 
and A 2 are 2|A||£/ A | 2 |F B | 2 and 2\X\\U A \ 2 \V B \ 2 respec- 
tively . All these products are shown to be small for 
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FIG. 2: a) Concurrence vs. dimensionless separation £ for 
r = txv/AQ. ~ A/8 and couplings K = K , WK , lOOifo 
and lOOOifo (bottom to top) b) Zoom around £ = 1 for the 
strongest coupling K = 1000-Kq- (Color online). 



the regions of interest discussed here, £ < 2. The same 
techniques can be extended to all orders in perturbation 
theory since the bounds to the different contributions can 
be grouped and treated as power series, giving rise to cor- 
rections that remain negligible as long as \A\, \X\, \Ua\ 2 
and \Vb\ 2 are small enough, like in the parameter range 
explored in this work. Finally, note that similar calcu- 
lations and results can be obtained in the case in which 
the qubits have close but different frequencies. 



III. ENTANGLEMENT DYNAMICS AND 
SINGLE PHOTONS 

We will use the concurrence C to compute the entan- 
glement of this state, which is given by 



C(pab) = - max{|p 2 3| 
c 



VPnPte , \pu\ - VP22P33 , 0} . 

(12) 

Since all quantities depend only on two dimensionless 
numbers, £ and K, we can perform a rather exhaustive 
study of the dynamics of entanglement between both 
qubits. To cover the widest possible spectrum of exper- 
iments, we have chosen coupling strengths over two or- 
ders of magnitude, K/K = 1, 10, 100, 1000. The smallest 
value K = 1.5 • 10~ 4 , which corresponds to g/n ~ 175 
MHz and fl/2ir ~ 10 GHz, that is for instance a charge 
qubit in the strong coupling limit with a transmission 
linens. The largest value, K = 1000 Kq corresponds to 
g ~ 2tt x 500 MHz and fl ~ 2tt x 2 GHz, and typically 
correspond s to a flux qubit directly coupled to a trans- 
mission lin d 19 ! 21 !, as shown in Fig. [lp-c. Note that in this 
case, by building the qubit directly on the line, much 
larger couplings can be achieved, in the ra nge of 1 — 4 
GHz, with 800 MHz being recently obtained^. 

In Fig. [2] we plot the value of the concurrence for two 
qubits which are separated a distance r = A/8, using the 
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FIG. 3: Concurrence (dash) and probability of excitation of 
atom B (solid) vs. dimensionless time, fit. Qubits are sep- 
arated by r — A/12 (circles) and A/8 (crosses) and have a 
coupling strength K = lOOOifo- Note that, following micro- 
causality, the excitation probabilities do not depend on the 
separation r outside the light cone (Color online). 



couplings discussed before. Note how the entanglement 
jumps discontinuously to a measurable value right inside 
the light cone (£ > 1), signalling the arrival of photons. 
Furthermore, even a certain amount of entanglement ap- 
pears outside the light cone, before photons could be ex- 
changed. This is best seen for the largest couplings, as 
Fig. [2]d illustrates. 

The dynamics looks even more exciting when we go 
back to lab time and space. Fig. [3] shows the concur- 
rence and the excitation probability of qubit B, ps = 
|^b| 2 / c + 0(d 4 ), for two different separations, r = A/12 
and r = A/8. The probability of excitation appears as 
independent of the qubit separation. This is exactly the 
case for the lowest order considered here, which only ac- 
counts for B self-interaction, and at all orders in per- 
turbation theory GO outside the light cone of this setup 
(region I in Fig. [1^,) . This is in full agreement with mi- 
crocausality. However, as can be seen in Fig. [3j what 
was a tiny concurrence jumps to a sizable value when 
crossing the light cone fit = 2n/12 and fit = 2n/8. In 
other words, from the experimental point of view, it is 
the entanglement between the qubits and not the excita- 
tion probability ps what best signals the presence of a 
light cone and a finite propagation speed. 



repeated many times, not only to gather sufficient statis- 
tics, but also to resolve different of instants of time be- 
fore and after the light-cone boundary. This may seem a 
daunting task, but thanks to the speed at which quantum 
circuits operate and their fast repetition rate, it will be as 
demanding as recent experiments realizing a controllcd- 
NOT gate^ or full two-qubit tomography^. 

The actual experimental challenge, though, arises from 
the need to perform quantum measurements of the qubit 
state and ensuring that this state is not altered by the 
ongoing dynamics. One possibility is to perform very fast 
measurements of the qubits, which means faster than 
I/O. The typical response of measurement apparatus, 
which in the case of SQUIDs is around a few nanoseconds, 
sets an upper limit on the qubit and photon frequencies 
of a few hundreds of megahertzs, though we expect this 
to be improved in the near future. 

Another more reliable approach is to connect and dis- 
connect the coupling between the qubit and the transmis- 
sion line. In this manner, we could prepare, entangle, and 
finally measure the qubits without interference or decay 
processes. If we work with flux qubits, a simple approach 
is to apply a very large magnetic flux on both qubits, 
taking the qubit away from its symmetry point. From 
a mathematical point of view, this amounts to adding a 
large contribution Ea x A B to the Hamiltonian. If done 
very quickly, the field projects the qubit on the same ba- 
sis on which the coupling operates, eliminating the pos- 
sibility of spontaneous emission. One would still need to 
combine the switching of this flux with short pulses that 
rotate the qubit basis in order to perform a complete set 
of measurements. The last and most elegant possibility 
is to effectively switch off all couplings between the qubit 
and the surrounding field. This can be achieved using a 
direct coupling between the qubit and the transmission 
line, with an scheme that incorporates an intermediate 
loop [Fig. []}:]. As we have shown in a recent work 2 ^, the 
result is a coupling that can be rotated and completely 
deactivated in a time of about 0.1 ns, that is the time 
needed to inject flux through the loop. The advantage is 
that, contrary to the case of a large external flux, the in- 
fluence of the line is completely suppressed and makes it 
possible to easily rotate the qubits to perform all needed 
measurements. 



CONCLUSIONS 



IV. EXPERIMENTAL IMPLEMENTATION 

In order to study the dynamics of quantum correla- 
tions between the two superconducting qubits, one has 
to perform a partial or full tomography of their state. In 
the first and simpler case, performing measurements in 
different basis should be enough to gather an entangle- 
ment quantifier, such as a Bell inequality violation or, as 
studied in this paper, the concurrence. This has to be 



Summing up, in this work we have proposed a circuit- 
QED experiment to study the dynamics of entanglement 
between two qubits that interact by exchanging traveling 
photons. Our work focuses on the existence of a finite 
propagation speed, the appearance of a light cone, the 
notion of microcausality and the possibility of achieving 
entanglement both by means of the correlated fluctua- 
tions of the vacuum and by photon exchange. The re- 
sulting predictions have a wide interest that goes beyond 
the assessment of microcausality in the QED of quan- 
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turn circuits, demonstrating that the open transmission 
line is a useful mediator of entanglement, much like cav- 
ities and zero-dimensional resonators. Furthermore, the 
experiment we propose is also among the simplest ones 
that can probe the effective QFT for waveguides, both as- 
serting the existence of propagating single photons and 
probing the dispersion relation at the single-photon level. 
Finally, we have shown that entanglement via traveling 
photons works better for stronger qubit-line couplings, 
making it one of the first potential applications of the 
ultrastrong coupling regime^. 
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Appendix 

In this appendix we will give further details on the 
computations of the relevant magnitudes \X\, \Ua\ 2 , 
\Vb\ 2 which are necessary to compute the concurrence 
Eq.Q. With Eqs. @-([7|) and the commutation rela- 
tions below Eq. ([I]): 



X = 



d 2 Nv 



ti 2 



dk\k\ {e lkr I t - 



^ — iki 



la) (A.l) 



with 



I na = ( dt 2 ,i f 2 ' 1 dti, 2 e in( * 2 - tl) e- lt,|fc|( * 2 ' 1 - tl - 2) (A.2) 
Jo Jo 

Notice that the term with I t i gives the non-RWA proba- 
bility amplitude associated to a single photon emission of 
qubit B followed by an absorption of qubit A. Performing 
the time integrations, inserting them in Eq. (A.l) and 



after some algebra, X can be given as a combination of 
integrals of the form: 



f 

Jo 

f 

Jo 



dk 



dk 



dk 



dk 



cos( £7) 
k + p 

cos( £7) 
k -/3 

sin( kj) 
k + p 

sin( k 7) 
k - P 



= — sin(7/3)s«(7/3) — cos(7/3)cz(7/3) 

= — sin(7/3)si(7/3) — cos(7/3)d(7/3) 

- 7rsin(7/3) (A.3) 

= sin(7/3)Cz(7/3) — cos(7/3)s«(7/3) 

= - sin(7/3)C?(7/3) + cos(7^)si(7/3) 

+ 7rcos(7/3) 



with 7, p > and the conventions irP^lfor the si and Ci. 
Putting all together we find: 



x = y(^p£sin(p) - eT((l 



zt_){C{t_)-S{t_) 
™(t_)6(1 - 0) + (1 - iT+)(-C(r + ) - S(t + ) 



-7rsin(r + ))+)(-r_ + i){-SC(r-) + CS(t-) + tt 
cos(r_)) + (-7+ - i)(-SC(T+) + CS(t+) +ttcos{t + )) 

-2) -e^((l-M-_)(-<7(T_) - S(r-) -7rsin(r_) 
9« - 1)) + (1 + ir+)(-C(r + ) - S(r+)) + (r_ + *) 
(SC(t_) - CS(t_)) + (t_ - i)(SC(T+) - CS(t+) + ir 
cos(t + )) - 2) - 2 - 2C(p) - 2S(p) - p(-2SC(p) 
+2CS( P )) (A.4) 



where £ has been defined in the main text, p = fir/t) 
is a dimensionless distance, t_ = p(l — £) = p — fit, 
t + = p(l + = p + fit, and we define C(x) = 
cos(x) Ci(x), S(x) = sin(a;) si(x), CS(x) = cos(x)si(x), 
SC(x) = sin(x) Ci(x). Notice the dependence with the 
spacetime region through the factors with the Heaviside 
function 0. 

Now we come to the emission probabilities |?7a| 2 , 
|Vb| 2 , which are given by 



\U A \ 2 = (Q\S+S^\0},\Vb\ 2 = (0\S b S+ 



(A.5) 



Following similar techniques we find that \Ua\ 2 = 
/+(nt), \V B \ 2 =/_(fit), where: 



f ± = y(7rf2i±2(cos(rU) + ntSi(flt) ~ 1)) (A.6) 



where Si must not be mistaken by si, Si = si + 7r/2 as 
usual. 



Finally, notice that A (Eq.Q) is a sum of two terms 
like I V a 1 2 aud \Vb\ 2 with the time ordering operator T 
and that Ua * Vb (Eq.(ll)) is similar to X without T. 
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